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Q . Exclusive semileptonic B decays to radially excited charmed mesons are 

*0 , investigated at the first order of the heavy quark expansion. The arising 

leading and subleading Isgur-Wise functions are calculated in the framework 
of the relativistic quark model. It is found that the 1/mq corrections play 
an important role and substantially modify results. An interesting interplay 



> 

l/"") ■ between different corrections is found. As a result the branching ratio for the 

B —* D'ev decay is essentially increased by 1/mq corrections, while the one 
for B — > D*'ev is only slightly influenced by them. 



On 

ON 

^ 1 I. INTRODUCTION 

a: 

Ph! The investigation of semileptonic decays of B mesons to excited charmed mesons rep- 

resents a task interesting both from the experimental and theoretical point of view. The 
current experimental data on the semileptonic B decays to the ground state D mesons in- 
dicate that a substantial part (~ 40%) of the inclusive semileptonic B decays should go to 
excited D meson states. First experimental data on some exclusive B decay channels to 
excited charmed mesons are becoming available now and more data are expected in 
near future. Thus the comprehensive theoretical study of these decays is necessary. The 
presence of the heavy quark in the initial and final meson states in these decays consider- 
ably simplifies their theoretical description. A good starting point for this analysis is the 
infinitely heavy quark limit, ttlq — > 00 [Q. In this limit the heavy quark symmetry arises, 
which strongly reduces the number of independent weak form factors |J. The heavy quark 
mass and spin then decouple and all meson properties are determined by light-quark degrees 
of freedom alone. This leads to a considerable reduction of the number of independent form 
factors which are necessary for the description of heavy-to-heavy semileptonic decays. For 
example, in this limit only one form factor is necessary for the semileptonic B decay to 
S- wave D mesons (both for the ground state and its radial excitations), while the decays to 
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P states require two form factors ||. It is important to note that the heavy quark symmetry 
requires that in the infinitely heavy quark limit matrix elements between a B meson and 
an excited D meson should vanish at the point of zero recoil of the final excited charmed 
meson in the rest frame of the B meson. In the case of B decays to radially excited charmed 
mesons this is the result of the orthogonality of radial parts of wave functions, while for the 
decays to orbital excitations this is the consequence of orthogonality of their angular parts. 
However, some of the 1/toq corrections to these decay matrix elements can give nonzero 
contributions at zero recoil. As a result the role of these corrections could be considerably 
enhanced, since the kinematical range for B decays to excited states is a rather small region 
around zero recoil. Recent calculations of semileptonic B decays to orbitally excited (P- 
wave) charmed mesons with the account of the 1/rriQ corrections support this observation 
|||7]]. Our calculations J7| in the framework of the relativistic quark model show that some 
rates of B decays to orbitally excited charmed mesons receive contributions from first order 
I/frig corrections approximately of the same value as a leading order contribution. In this 
paper we extend our analysis to B decays to radially excited D mesons. 

Our relativistic quark model is based on the quasipotential approach in quantum field 
theory with a specific choice of the quark-antiquark interaction potential. It provides a 
consistent scheme for the calculation of all relativistic corrections at a given v 2 /c 2 order and 
allows for the heavy quark 1/mg expansion. In preceding papers we applied this model to 
the calculation of the mass spectra of orbitally and radially excited states of heavy-light 
mesons ||, as well as to a description of weak decays of B mesons to ground state heavy 
and light mesons PJlCj. The heavy quark expansion for the ground state heavy-to-heavy 



semileptonic transitions [|TTJ has been found to be in agreement with model-independent 
predictions of the HQET. 

The paper is organized as follows. In Sec. II we cary out the heavy quark expansion for the 
weak decay matrix elements between a B meson and radially excited charmed meson states 
up to the first order in 1/toq using heavy quark effective theory (HQET). In our analysis we 
follow HQET derivations for the matrix elements between ground states [|12| , P~3f and for the 
matrix element between a B meson and orbitally excited charmed meson ||, as well as a 
general analysis of these matrix elements in Ref. |TJ]]. In Sec. Ill we describe our relativistic 
quark model. The heavy quark expansion for decay matrix elements is then carried out up 
to the first order I/uiq corrections and compared to the model-independent HQET results 
in Sees. IV- V. We determine the leading and subleading Isgur-Wise functions and give our 
predictions for decay branching ratios in the heavy quark limit and with the account of 1 / tuq 
corrections. The electron spectra for the considered decays are also presented. Section IV 
contains our conclusions. 



II. DECAY MATRIX ELEMENTS AND THE HEAVY QUARK EXPANSION 

The matrix elements of the vector current (JY = cry^fr) and axial vector current (J^ = 
C7/x75&) between B and radially excited D' or D*' mesons can be parameterized by six 
hadronic form factors: 

^m D ,m B 
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(D'{v')\cYH>\B{v)) = 
(D*'(v', e)\cq^b\B{v)) 



y/m D *im B 
(D*'{v',e)\c^b\B(v)} 
y/m B */m B 



h Al (w + l)e*" - (h A y + h M v'»){e* ■ v), 



(1) 



where v (V) is the four- velocity of the B (£)(*)') meson, is a polarization vector of the final 
vector charmed meson, and the form factors hi are dimensionless functions of the product 
of velocities w = v ■ v '. The double differential decay rates expressed in terms of the form 
factors are 



d 2 r D , 

dwdcos 9 

d 2 ry, 

dwdcos 9 



3T r 3 (w 2 - 1) 3/2 sin 2 9 [(1 + r)h+ - (1 - r)h 



sin 2 6* (w — r*)h Al + (w 2 — l){h A3 + r*h A2 ) 



+ (l - 2r*w + r 



*2^ 



1 + cos 2 9)[h 2 Al + (w 2 - l)h 2 v }- 4 cos 9Vw 2 - lh Al h v k (2) 



where T = G 2 F \V c b\ 2 m 5 B / (192n 3 ) , r = mD>/m B , r* = m D *i/m B and 9 is the angle between 
the charged lepton and the charmed meson in the rest frame of the virtual W boson. 

Now we expand the form factors hi in powers of 1/rriQ up to first order and relate 
the coefficients in this expansion to universal Isgur-Wise functions. This is achieved by 
evaluating the matrix elements of the effective current operators arising from the HQET 
expansion of the weak currents. For simplicity we limit our analysis to the leading order in 
a s and use the trace formalism ||15|| . Following Ref. [H], we introduce the matrix 



i+ i 



p7 7m - p vl5 



(3) 

| doublet [] with 

four-velocity v. At leading order of the heavy quark expansion (mq — > oo) the matrix 
elements of the weak current between the ground and radially excited states destroyed by 
the fields in H v and H' v , respectively, are given by 



composed from the fields P v and P*^ that destroy mesons in the j J 



en - t$rh® = £ (n) MTr{^:,r#„}, 



(4) 



where hffi' is the heavy quark field in the effective theory. The leading order Isgur-Wise 
function ^ n \w) vanishes at the zero recoil (w = 1) of the final meson for any T, because of 
the heavy quark symmetry and the orthogonality of the radially excited state wave function 
with respect to the ground state one. 

At first order of the 1 / uiq expansion there are contributions from the corrections to the 
HQET Lagrangian 



1 Here j is the total light quark angular momentum, and the superscript P denotes the meson 
parity. 
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5C 



2m, 



AQ) - 

'-'1m — 



2m, 



(5) 



and from the tree-level matching of the weak current operator onto effective theory which 
contains a covariant derivative D x = d x — ig s t a A x 



cYb ^ h^) (v - ^-$Y + 1 



2m c 2mfe 

The matrix elements of the latter operators can be parameterized as 



(6) 



h<$iD x Yh® = Tr{^H v ,YH v }, 
h$Til} x hf = Tr{^ ] H v ,YH v }. 



The most general form for £^ is [|12 



(Q) 



(7) 



(8) 



The equation of motion for the heavy quark, i{v ■ D)h^ = 0, yields the relations between 
the form factors £^ 



(<0, 



.(c) 
,3 
: (b) 




0. 



(9) 



The additional relations can be obtained from the momentum conservation and the defi- 
nition of the heavy quark fields which lead to the equation id„(h$Y h^) = (Av v — 
AWt/J^r^, implying that 



.(c) 



.(b) 



er + cr = - A (n M)e (n) 



(10) 



Here A(A^ n ^) = M(M^) — vtlq is the difference between the heavy ground state (radially 
excited) meson and heavy quark masses in the limit ijiq — ► oo. This equation results in the 
following relations 



(b) 



0. 



(11) 



The relations (0) and ( ]TT|) can be used to express the functions in terms of |s(= £3^ 
—£3^) and the leading order function £( n ); 



£ 



(c) 



.(&) 



'AW lAW-A 

+ 2 w - 1 
'A lA^-A' 
2" ~ 2 w - 1 
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Jc) 



AM iaW + A v 



(6) 



2 ' 2 w + 1 
'A lA^+AN 



w + 1 
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2 2 w + 1 



(n) 



1 



w + 1 



(12) 



The matrix elements of the I/tuq corrections resulting from insertions of higher- 
dimension operators of the HQET Lagrangian ([5]) have the structure [[R] 

dxT{£<&(x) [J$Th®] (0)} = 2 X [ c) Tr {h v ,TH v } + 2Tr ^H^^-TH^ , 

dxT{C?l(x) [htfrhg*] (0)} = 2 X S 6) Tr {S V ,TH V } + 2Tr L%H v ,r^ia a P H v } . (13) 

The corrections coming from the kinetic energy term Okm do not violate spin symmetry and, 
hence, the corresponding functions Xi effectively correct the leading order function 
The chromomagnetic operator O mag , on the other hand, explicitly violates spin symmetry. 
The most general decomposition of the tensor form factor XaB ^ s E^jH 



(c) (c) (c) . 

Xa/3 = X2 Val/3 ~ X 3 W a/3, 
(6) (b) , (b). 



(14) 



The functions xt contribute to the decay form factors (|]) only in the linear combination 

Xb = 2x? — 4(w — 1)X2 + 12x3^. Thus five independent functions £ 3 , Xb and x«(= Xi )j as 
well as two mass parameters A and A^ n ' are necessary to describe first order 1/mg corrections 
to matrix elements of B meson decays to radially excited D meson states. The resulting 
structure of the decay form factors is 



h. 



+ s c [2xi - 4(w - 1)X2 + 12x 3 ] + e bXb , 



2| 3 - A (n) + 



W — 1 



Sb 



2& - A 



AW - A s 

W — 1 



v = e (n) + ^ 



2X1-4X3+ A< n > + 



A (n) - A ' 
w — 1 



+£b 



X6+ A 



AM - A s 



hA 2 



+£b |x& + 



2xi - 4X3 + 
w-1 , A 



w — 1 

w — 1 



e (n) - 2f 3 



w + 1 
2 

w + 1 



W + 1 

AH - A 



A (n) + 



AW - A s 
w — 1 



A (n) + 



W — 1 / 

AW - A 



£(») _ 2 | 3 



2xi - 4X2 - 4x3 + 



w — 1 
w — 1 



£ {n) + 6 



w + 1 



A (n) + 



AW - A s 
w — 1 



w + 1 
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(15) 



where e Q = l/(2m Q ). 

In the following sections we apply the relativistic quark model to the calculation of 
leading and subleading Isgur-Wise functions. 



III. RELATIVISTIC QUARK MODEL 

In the quasipotential approach, a meson is described by the wave function of the bound 



quark- ant iquark state, which satisfies the quasipotential equation [16fl of the Schrodinger 
type [0: 



'b 2 {M) 



*Af(p) 



where the relativistic reduced mass is 



d 3 q 



y(p,q; M)tf M (q) 



(16) 



M 4 — 



[m\ - m 2 Q ) 2 



4M 3 



(17) 



Here m q ^Q are the masses of light and heavy quarks, and p is their relative momentum. In 
the center-of-mass system the relative momentum squared on mass shell reads 



b 2 {M) 



[M 2 - (m q + m Q ) 2 } [M 2 - (m q - m Q f 
4M 2 



(18) 



The kernel V(p, q; M) in Eq. ( Pp is the quasipotential operator of the quark- ant iquark 
interaction. It is constructed with the help of the off-mass-shell scattering amplitude, pro- 
jected onto the positive energy states. An important role in this construction is played by 
the Lorentz-structure of the confining quark- ant iquark interaction in the meson. In con- 
structing the quasipotential of the quark-antiquark interaction we have assumed that the 
effective interaction is the sum of the usual one-gluon exchange term and the mixture of 
vector and scalar linear confining potentials. The quasipotential is then defined by fi~8 



V(p, q; M) = u q (p)u Q {-p)V{p, q; M)u q {q)u Q {-q) 

[4 

= u q {p)u Q (-p)\ -a s LV( k )7g7Q 



+V^MT q T Q ,, + KonfO) \u g (q)u Q (-q), 



(19) 



where a s is the QCD coupling constant, is the gluon propagator in the Coulomb gauge 
and k = p — q; 7 M and u(p) are the Dirac matrices and spinors 



u x (p) 



\ 



e(p) + m | 




Hp) ' 


crp J 


\ e(p)+m / 



X 



(20) 
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with e(p) = Vp 2 + m 2 - The effective long-range vector vertex is given by 



r M (k) = 7 M + <^ ( V r > 



(21) 



where k is the Pauli interaction constant characterizing the nonperturbative anomalous chro- 
momagnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 



Ko» = (i 



■ £ 



:) (Ar + B), 



reproducing 



V^(r) = e(Ar + B) 



Konf (r) = Konf(0 + Konf W = ^ + 5, 



(22) 



(23) 



where e is the mixing coefficient. 

The quasipotential for the heavy quarkonia, expanded in v 2 /c 2 , can be found in 
Refs. |T8|,19| and for heavy-light mesons in ||. All the parameters of our model, such as 
quark masses, parameters of the linear confining potential, mixing coefficient e and anoma- 
lous chromomagnetic quark moment k, were fixed from the analysis of heavy quarkonia 
18[| and radiative decays ||20|| . The quark masses rrib = 4.88 GeV, m c = 1.55 GeV, 
0.33 GeV and the parameters of the linear potential A = 0.18 GeV 2 



masses 
m s = 0.50 GeV, m u a 



and B = —0.30 GeV have the usual quark model values. In Ref. j! 1] we have considered 
the expansion of the matrix elements of weak heavy quark currents between pseudoscalar 
and vector meson ground states up to the second order in inverse powers of the heavy quark 
masses. It has been found that the general structure of the leading, first, and second or- 
der 1/mq corrections in our relativistic model is in accord with the predictions of HQET. 
The heavy quark symmetry and QCD impose rigid constraints on the parameters of the 



long-range potential in our model. The analysis of the first order corrections [III] allowed 
us to fix the value of the Pauli interaction constant k — — 1. The same value of k was 
found previously from the fine splitting of heavy quarkonia 3 Pj- states | JL8|| . Q Note that 
the long-range chromomagnetic spin-dependent interaction in our model is proportional to 
(1 + k) and thus vanishes for n = — 1 in agreement with the flux tube model p2| . The value 
of the mixing parameter of vector and scalar confining potentials e — — 1 has been found 
from the analysis of the second order corrections [|Tl |. This value is very close to the one 
determined from radiative decays of heavy quarkonia |^0| . 



In order to calculate the exclusive semileptonic decay rate of the B meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current J w = C7^(l — 7 5 )6 
between a B meson and a radially excited D^' meson takes the form 



(DW\J?(0)\B) 



d 3 p d 3 q 
(2vr) 6 



*£»w(p)r«(p,q)*B(q), 



(24) 



2 It has been known for a long time that the correct reproduction of the spin-dependent part of 
the quark-antiquark interaction requires either assuming the scalar confinement or equivalently 
introducing the Pauli interaction with k = —1 [21.18,19] in the vector confinement. 
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where r M (p, q) is the two-particle vertex function and ^ B D {*), are the meson wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
frame. The contributions to T come from Figs. 1 and 2.[] In the heavy quark limit m& jC — > oo 
only T^ 1 ' contributes, while r^ 2 ) contributes at 1/rriQ order. They look like 

r«(p, q) = u c ( Pc )l,(l - 7 5 )^(g6)(27r) 3 5( P(? - q q ), (25) 

and 

F^(p, q) = U c (p c )u q (p q ){ lQtl (l - ^)_±Xj_^ V (p fl - q 9 ) 

+V(P ? - - 7S)}«6(ft)«,(«,), (26) 

where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, Q = c or b, k = 
p c - A; k' = q 6 + A; A = p DW , - p B , 

A (_ } = e(p) - (m 7 + 7°(7p)) 
2e(p) 

Here @ 

3 



3 



and are three four-vectors given by 



n W"(p) = j— , <k- + r 1 , £ = Vp 2 + M 2 . 

UJ [M' 3 M{E + M)y VP 

It is important to note that the wave functions entering the weak current matrix element 
(^4f) are not in the rest frame in general. For example, in the B meson rest frame, the 
D^' meson is moving with the recoil momentum A. The wave function of the moving D^' 
meson \l/ £>(.)' a is connected with the D^' wave function in the rest frame $ fl (»)/ by the 



transformation p3 



«Wa(p) = £ c V WJ^WJ*d(*Vo(p), (27) 

where R w is the Wigner rotation, La is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D 1 ^ 2 (R) in spinor representation is given by 



3 The contribution I^ 2 ) is the consequence of the projection onto the positive-energy states. Note 
that the form of the relativistic corrections resulting from the vertex function I^ 2 ) is explicitly 
dependent on the Lorentz structure of the gg-inter action. 



S 



(J °)<WJ = S- 1 (p c>q )S(A)S(p), (28) 

where 

^ V 2m ( e(p) + rn) 
is the usual Lorentz transformation matrix of the four-spinor. 



IV. LEADING AND SUBLEADING ISGUR-WISE FUNCTIONS 

Now we can perform the heavy quark expansion for the matrix elements of B decays 
to radially excited D mesons in the framework of our model and determine leading and 
subleading Isgur-Wise functions. We substitute the vertex functions and given by 
Eqs. (pB[) and fl2~6|) in the decay matrix element (J2||) and take into account the wave function 
transformation (^7]). The resulting structure of this matrix element is rather complicated, 
because it is necessary to integrate both over d 3 p and d 3 q. The 5 function in expression 
( pip permits us to perform one of these integrations and thus this contribution can be easily 
calculated. The calculation of the vertex function r^ 2 ^ contribution is more difficult. Here, 
instead of a 5 function, we have a complicated structure, containing the Qq interaction 
potential in the meson. However, we can expand this contribution in inverse powers of 
heavy (b, c) quark masses and then use the quasipotential equation in order to perform one 
of the integrations in the current matrix element. We carry out the heavy quark expansion 
up to first order in I/tuq. It is easy to see that the vertex function contributes already 
at the subleading order of the 1 / txiq expansion. Then we compare the arising decay matrix 
elements with the form factor decomposition ([J) and determine the corresponding form 
factors. We find that, for the chosen values of our model parameters (the mixing coefficient 
of vector and scalar confining potential e — — 1 and the Pauli constant k — — 1), the resulting 
structure at leading and subleading order in I/tuq coincides with the model-independent 
predictions of HQET given by Eq. ([15]). We get the following expressions for leading and 
subleading Isgur-Wise functions: 

:<!),.„•, - fJ_V f 1V T(0, 2J, A ,(„), 



«.<«) = (— m. + 5— ) (l + 3^Tl) «*° «■ (30) 

—l) J (2^^H P+ M p ,.,,( tt , + l) A j^ (p) ' (31) 

a -ssTT^r^w- (32) 

80 w + 1 w — 1 
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AW / 2 y/ 2 r d*p - {l)sd ( 2e q \ (0) 



2 f d 3 p - (0 ) / , 2e 9 



»(«0 = A ( -- ) / t^«v( P + U "7. , Tl A ] # 4 (P) - 3^ 1)sd (p) . (34) 



.10 + 1/ 7 (27r) 3ri)( * , \^ M D (.),(w+l 

Here we used the expansion for the S'-wave meson wave function 

^ M = + ^mS Q (i^T + d M ^S Sd ) + 0(l/m|), 

where ip^ is the wave function in the limit mq — > oo, ' and ip$ sd are the spin- 
independent and spin-dependent first order 1/mg corrections, dp = —3 for pseudoscalar 
and dy = 1 for vector mesons. The symbol = in the expressions (j3l"l)-(|3~4f) for the sublead- 
ing functions Xi{w) means that the corrections suppressed by an additional power of the 
ratio (w — l)/(w + 1), which is equal to zero at w = 1 and less than 1/6 at u> ma x, were 
neglected. Since the main contribution to the decay rate comes from the values of form 
factors close to w = 1, these corrections turn out to be unimportant. 

It is clear from the expression (|29D that the leading order contribution vanishes at the 
point of zero recoil (A = 0, w — 1) of the final D^' meson, since the radial parts of the 
wave functions ^ £>(*)' and ^ b are orthogonal in the infinitely heavy quark limit. The 1/tuq 
corrections to the current ( |T2"D also do not contribute at this kinematical point for the same 
reason. The only nonzero contributions at w = 1 come from corrections to the Lagrangian [] 
X3( w ) an d Xbiw). From Eqs. flT5| ) one can find for the form factors contributing to 
the decay matrix elements at zero recoil 

Ml) = [2xi(l) + 12X3(1)] + ^6X6(1), 

h Al (l) = e c [2xi(l) - 4x 3 (l)] + £6X6(1)- (35) 

Such nonvanishing contributions at zero recoil result from the first order 1/mq corrections 
to the wave functions (see Eq. Q3"4] ) and the last terms in Eqs. fl3"T|), (|3"3"D). Since the kine- 
matically allowed range for these decays is not broad ( 1 < w < u> max ~ 1.27) the relative 
contribution to the decay rate of such small 1/rriQ corrections is substantially increased. 
Note that the terms £q(Aw — A)^ n \w) / (w — 1) have the same behaviour near w = 1 as 
the leading order contribution, in contrast to decays to the ground state D^*' mesons, where 
1/rriQ corrections are suppressed with respect to the leading order contribution by the fac- 
tor (w — 1) near this point (this result is known as Luke's theorem |12|]). Since inclusion of 
first order heavy quark corrections to B decays to the ground state D^*' mesons results in 
approximately a 10-20% increase of decay rates |Tl|Jl3[| , one could expect that the influence 



of these corrections on decay rates to radially excited D^*>' mesons will be more essential. 
Our numerical analysis supports these observations. 



There are no normalization conditions for these corrections contrary to the decay to the ground 
state Z)W mesons, where the conservation of vector current requires their vanishing at zero recoil 
& 
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V. NUMERICAL RESULTS AND PREDICTIONS 



In Table | we present the masses of radially excited D' and D*' as well as mass parameters 
A calculated in the framework of our model ||. Our prediction for the D*' mass is in good 
agreement with the DELPHI measurement |24|] . Other radially excited states have not been 



observed yet. Thus we use our predictions for numerical calculations. The values of leading 
and subleading Isgur-Wise functions (|29|)-(|34D and their slopes at the point of zero recoil of 
the final meson are given in Table [TT1 In Fig. 3 we plot our results for the leading order 
Isgur-Wise function ^'(w) and the current correction function £3(10). The functions Xi(w), 
X2(w), X3(w), and Xb(w) are plotted in Fig. 4. We see that the functions, parameterizing 
chromomagnetic corrections to the HQET Lagrangian, are rather small in accord with the 
HQET based expectations. 

We can now calculate the decay branching ratios by integrating double differential decay 
rates in Eq. (D). Our results for decay rates both in the infinitely heavy quark limit and 
taking account of the first order 1/mq corrections as well as their ratio 

Br(B ^ D^'ev) withl/mQ 

are presented in Table TT|. We find that both 1/rriQ corrections to decay rates arising 
from corrections to HQET Lagrangian ([H])-(34), which do not vanish at zero recoil, and 



corrections to the current (pOf), fll2"|) , vanishing at zero recoil, give significant contributions. 
In the case of B — > D'ev decay both types of these corrections tend to increase the decay 
rate leading to approximately a 75% increase of the B — > D'ev decay rate. On the other 
hand, these corrections give opposite contributions to the B — > D*'ev decay rate: the 
corrections to the current give a negative contribution, while corrections to the Lagrangian 
give a positive one of approximately the same value. This interplay of 1/ mg corrections only 
slightly (~ 10%) increases the decay rate with respect to the infinitely heavy quark limit. 
As a result the branching ratio for B — > D'ev decay exceeds the one for B — > D*'ev after 
inclusion of first order I/tuq corrections. In the infinitely heavy quark mass limit we have 
for the ratio Br(B — > D'ev)/Br(B — > D*'ev) = 0.75, while the account of 1/rriQ corrections 
results in the considerable increase of this ratio to 1.22. 

In Table [TTI| we also present the sum of the branching ratios over first radially excited 
states. Inclusion of 1/ uiq corrections results in approximately a 40% increase of this sum. 
We see that our model predicts that 0.40% of B meson decays go to the first radially excited 
D meson states. If we add this value to our prediction for B decays to the first orbitally 
excited states 1.45% @, we get the value of 1.85%. This result means that approximately 
2% of B decays should go to higher excitations. 

In Figs. 5 and 6 we plot the electron spectra (l/r )(dr/d?/) for B — > D'ev and B — > D*'ev 
decays. Here y = 2E e /rriB is the rescaled lepton energy. These differential decay rates can 
be easily obtained from double differential decay rates (||), using the relation y = 1 — rw — 
r\Jw 2 — 1 cos# and then integrating in w over [(1 — y) 2 + r 2 ]/[2r(l — y)] < w < (1 +r 2 )/ (2r). 
We present our results both in the heavy quark limit rriQ — > 00 (dashed curves) and with the 
inclusion of first order 1/mg corrections (solid curves). From Fig. 6 we see that inclusion of 
1/tuq corrections significantly changes the shape of electron spectrum for B — > D*'ev decay. 
The maximum is considerably shifted to higher lepton energies. 
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VI. CONCLUSIONS 



In this paper we have carried out the heavy quark expansion for the decay matrix elements 
of weak currents between the B meson and radially excited £)W' meson states up to first 
order. It is found that five additional functions of the product of velocities w are necessary 
to parameterize first order I/tuq corrections. One of these functions ^(w) arises from 
corrections to the weak current. The other four functions Xi{ w ) (* = 1 • • • 3, 6) parameterize 
corrections to the HQET Lagrangian. The leading order function vanishes at the point 
of zero recoil (w = 1) of the final Z)W' meson due to the heavy quark symmetry and 
orthogonality of the radial parts of meson wave functions in the heavy quark limit. The 
contributions to the decay matrix elements coming from the corrections to the current also 
vanish at zero recoil for the same reason. Thus the only nonzero contributions to the weak 
decay matrix elements at w — 1 come from the corrections to the Lagrangian Xi,3,6 ( see 
Eq. Q55D), since there is no condition requiring them to vanish at this point as in the case 
of B decays to ground state D mesons. 

Then we apply the relativistic quark model for the consideration of semileptonic B — ► 
D^'eu decays. It is found that our model correctly reproduces the structure of decay ma- 
trix elements found from the heavy quark symmetry analysis. This allows us to determine 
the leading and subleading Isgur-Wise functions for this transitions. We find that both the 
relativistic transformation of the meson wave function from the rest frame to the moving 
one as well as the first order 1/rriQ corrections to meson wave functions give essential con- 
tributions to the subleading order functions. Thus, the account for corrections to the wave 
functions gives contributions to decay matrix elements which do not vanish at zero recoil. 
These contributions turn out to be rather small numerically. However, their role is consid- 
erably increased since the kinematical range for B decays to radially excited D^' mesons is 
small and thus the leading order contribution, vanishing at zero recoil, is suppressed. An- 
other important contribution to decay rates comes from the terms Eq(A^ — A)^ 1 ^ /(w — 1) 
originating from the corrections to the current. These terms turn out to be numerically 
important. We find an interesting interplay of these two types of I/tuq corrections. They 
contribute to the B — > D'ev decay rate with the same sign, but their contributions to the 
B — > D*'ev rate have opposite signs. As a result the former decay rate is substantially 
(1.75 times) increased by the inclusion of first order I/trq corrections while there is only a 
slight (1.1 times) increase of the latter decay rate. This leads to the increase in the ratio 
Br(B — > D'ev)/Br{B — > D*'ev) from 0.75 in the heavy quark limit to 1.22, when first order 
I/tuq corrections are taken into account. Finally, we find that the semileptonic B decays to 
first radial excitations of D mesons acquire in total 0.4% of the B decay rate. 
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TABLES 

TABLE I. Masses of radially excited £)(*)' mesons and the mass parameters A in our model. 



Parameter 


State 


Value (GeV) | 


m D i 


£>'(25 ) 


2.579 


m D *i 


D*'(2Si) 


2.629 


A 


B,D(1S) 


0.51 


Ad) 


D(2S) 


0.94 


m D / 


D' S (2S ) 


2.670 


m Dr 


D* S '(2S 1 ) 


2.716 


A s 


B S ,D S (1S) 


0.61 


A« 


D S (2S) 


1.05 



Exp. (GeV) |4 



2.637(9) 



TABLE II. Leading and subleding Isgur-Wise functions and their slopes p^. = — Jt^j^ 



i a 



and ^ = -5&K5Xi 
and defined ^^(w) 
(AW + A)/2. 



u)=l 



at zero recoil. We factored out (u; — 1) from the leading order form factor 
(w — l)S(u;). The values of the functions £3(1), Xi(l) are given in units 



E(w) 



Xi(w) 



X2{W) 



X3{W) 



Xb{w) 



Value at w = 1 2.2 
Slope at w = 1 2.6 



0.21 
-3.3 



0.18 
1.9 



-0.054 
3.1 



-0.023 
1.0 



-0.098 
2.1 



TABLE III. Decay rates T (in units of |V r c6 /0.04| 2 x 10" 15 GeV) and branching ratios BR (in 
%) for B (B s ) decays to radially excited D^' (£>i* }/ ) mesons in the infinitely heavy quark mass 
limit and taking account of first order 1/mg corrections. T,(B — > D^'eu) and T,(B S — ► D^'ev) 
represent the sum over the channels. R is a ratio of branching ratios taking account of 1/toq 
corrections to branching ratios in the infinitely heavy quark mass limit. 









-> OO 




With 1/toq 






Decay 


r 


Br 


L 


Br 


R 


B 


-> D'ev 


0.53 


0.12 


0.92 


0.22 


1.74 


B 


-> £>*'ez/ 


0.70 


0.17 


0.78 


0.18 


1.11 


E(B 


-» D^'ev) 


1.23 


0.29 


1.70 


0.40 


1.37 


B s 


-» D' s ev 


0.66 


0.16 


1.18 


0.28 


1.80 


B s 


-> -D s *'ei^ 


0.86 


0.20 


0.95 


0.22 


1.10 


E(B S 


- z^'H 


1.52 


0.36 


2.13 


0.50 


1.40 
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FIGURES 




FIG. 1. Lowest order vertex function contributing to the current matrix element (25). 
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FIG. 2. Vertex function r^ 2 ^ taking the quark interaction into account. Dashed lines correspond 
to the effective potential V in ( |I9D . Bold lines denote the negative-energy part of the quark 
propagator. 
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1 1.1 1.2 1.3 1.4 



FIG. 3. Isgur-Wise functions £W(w) (solid curve) and £3(10) (bold curve, in units (A^ 1 ) +A)/2) 
for the B — > D^'eu decay. 



16 




1 1.1 1.2 1.3 1.4 



FIG. 4. Isgur-Wise functions Xi( w ) (solid curve), X2(w) (bold curve), X3( w ) (dotted curve), 
and Xb{w) (dashed curve) for the B —> D^*>'eu decay in units (AW + A)/2. 
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FIG. 5. Electron spectra (1/To) (dT/dy) for the B — > D'eu decay as a function of the rescaled 
lepton energy y = 2E e /rriB- Dashed curves show the tuq — > oo limit, solid curves include first 
order 1/mg corrections. 
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FIG. 6. Electron spectra (1/To) (aT/dy) for the B — > D*'ev decay as a function of the rescaled 
lepton energy y = 2E e /rriB- Dashed curves show the tuq — > oo limit, solid curves include first 
order 1/toq corrections. 
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